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We give a stochastic extension of the Brane Calculus, along the lines of recent work by Cardelli and 
Mardare (8). In this presentation, the semantics of a Brane process is a measure of the stochastic 
distribution of possible derivations. To this end, we first introduce a labelled transition system for 
Brane Calculus, proving its adequacy w.r.t. the usual reduction semantics. Then, brane systems 
are presented as Markov processes over the measurable space generated by terms up-to syntactic 
congruence, and where the measures are indexed by the actions of this new LTS. Finally, we provide 
a SOS presentation of this stochastic semantics, which is compositional and syntax-driven. 

1 Introduction 

The Brane Calculus (BC) is a calculus of mobile processes designed for modeling membrane in- 
teractions within a cell. A process of this calculus represents a system of nested membranes, carrying 
their active components on membranes, not inside them. Membranes interact according to three reaction 
rules, corresponding to phagocytosis, endo/exocytosis, and pinocytosis. 

In the original definition, reaction rules do not consider quantitative aspects like rates, volumes, 
etc. However, it is important to address these aspects, e.g. for implementing stochastic simulations, or 
for connecting Brane Calculus with quantitative models at lower abstraction levels (such as stochastic 
TT-calculus and fc-calculus for protein interactions). 

In this paper, we introduce a stochastic semantics for the Brane Calculus. Clearly, a stochastic 
calculus could be obtained just by adding rates to reaction rules; however, the resulting "pointwise" 
rated reduction semantics is not fully satisfactory for several reasons. First, it is not compositional, i.e., 
reaction rates of a process are not given in terms of the rates of its components. Secondly, stochastic 
reaction rules are not easy to deal with in presence of large populations of agents (as it is often the case 
in biological systems), because we have to count large number of occurrences for calculating the effective 
reaction rates. Third, it does not generalize easily to other quantitative (e.g. geometric) aspects. 

To overcome these issues, we adopt a novel approach recently introduced by Cardelli and Mardare 
[ 8 ] , which is particularly suited when a measure of similarity of behaviours is important (similar ideas 
have been proposed for probabilistic automata ifTOl [T4l . and Markov processes JU 13 [13). The main 
point of this approach is that the semantics of a process is a measure of the stochastic distribution of 
the possible outcomes. Thus, processes form a measurable space, and each process is given an action- 
indexed family of measures on this space. For an action a, the measure \i a associated to a process P 
specifies for each measurable set S of processes, the rate jU a (5) 6 M + of a-transitions from P to (elements 
of) S. The resulting structures, called Markov processes (MPs), are not continuous-time Markov chains 
because each transition is from a state to a possibly infinite class of states (closed to the congruence 
relation over processes) and consequently cannot be described in a pointwise style. An advantage of 
this approach is that we can apply results from measure theory for solving otherwise difficult issues, 
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like instance-counting problems; moreover, process measures are defined compositionally, and can be 
characterized also by means of operational semantics in GSOS form. Finally, other measurable aspects 
of processes (e.g., volumes) can be dealt with along the same lines. 

However, the approach of (H has been applied to CCS only, and in order to adapt it to BC, we 
have to solve some problems. First, in order to define systems of BC as Markov processes, we need 
to define actions for brane systems, which corresponds to define a labelled transition system (LTS) for 
BC. Defining correctly a labelled transition system for a complex calculus like BC is notoriously a 
difficult task, because labels should describe precisely how a system can interact with the surrounding 
environment. In our case, we have been inspired by the so-called IPO construction ifTTl . following 
the approach of |[T3l l3ll for the case of Mobile Ambients. As a result, in Section [3] we introduce the first 
labelled transition system for BC, and prove its adequacy with respect to the original reduction semantics. 

A peculiar feature of this LTS is that after a transition, a system can yield a higher-order term, i.e., 
a system with "holes" (like 7T-calculus' "abstractions"). This has several consequences on our work. 
First, we have to define a suitable syntax for these higher-order terms; for this reason, in Section [2] we 
introduce a simply typed version of brane calculus from the beginning, extended with metavariables, X- 
abstractions and applications. Well-formedness of terms is guaranteed by a suitable typing system; e.g., 
membranes and systems are represented by ground terms of type mem and sys, respectively. Secondly, 
the bisimulation definition has to be accommodated in order to deal with transitions yielding higher-order 
terms (Section[3]>. Finally, also the definition of Markov kernel is affected, as we will see in Section|4} the 
kernel cannot be defined simply on the space of brane systems, but must consider also these higher-order 
terms — in fact, the measure will be defined over all well-typed terms, also those higher-order. 

After that these issues have been addressed and the Markov kernel has been defined, we can look for 
a simpler presentation of the semantics of Markov processes. In Section [5] we present a SOS system for 
processes, capturing the Markov kernel over processes: the stochastic bisimilarity induced by this SOS 
semantics corresponds to the Markov bisimilarity defined in Section [4] Therefore, this semantics can be 
fruitfully used for simulations, or for verifying system equivalences. 

Some concluding remarks and directions for further work are in Section [6j 

2 Brane Calculus 

In this section we recall Cardelli's Brane Calculus [7] focusing on its basic and finite version (without 
communication primitives, molecular complexes and replication). Here we adopt an alternative presenta- 
tion of the calculus: instead of the traditional term grammar of Q, we introduce an unstructured version 
but equipped with a type system. Although typed terms may seem unnecessary now, they will be useful 
in Section|3j where a labelled transition semantics will be introduced. 

Syntax The grammars for terms and types are specified below. 

{Terms) M::=0 | « | X \ a.M \ M\M \ M°M | MflMD 

«::=»„ | | S„ | s£ | ®„{M) 

{Types) t::=sys | mem | act 

The subscripted names n are taken from a countable set A, while term variables X are taken from a 
countable set 3£, assumed to be disjoint from each other. 
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F(X)=t a G {»„,»„, Q^} ae{^,@„} rhM:mem 

(var) — (act) — ; : (act-are) 

ThX:t Tha:act Tha(M):act 

n h a : act r 2 hM:mem rihM:mem r 2 h./V:mem 

(zero) — — ~ — : (a-pref) — — — — : ; (par) 



— (zero) — — — — : (a-pret) — — — — ; ; 

rhO:mem r 1 ,r 2 ha.M:mem ri,r 2 \~M\N : mem 

H h M : mem r 2 h N : sys r, h M : sys r 2 h N : 

r h o : sys ( ™ ld) ri,r 2 r-Md^VD :sys (1 ° C) ri,r 2 hMotf : sys 

Table 1 : Typing system for Brane Calculus 



Types are assigned to terms as usual. A type environment Y is a finite map from term variables to 
types. If r is an environment and X a variable not in the domain of T, we denote by T,X : t the environ- 
ment which assigns type T(Y) to each variable Y G dom{T) and type t to X; if Ti and T 2 have disjoint 
domains, Ti , T 2 denotes the environment which assigns type F\ (X) to variables X G dom{T\ ) and T 2 (F) 
to variables Y G Jom(r 2 ). The type inference rules are given in Table [T] Notice that this type system ad- 
mits only linear terms, that is, each variable can occur at most once. Indeed, in rules (a-pref), (par), (loc), 
and (comp), the environment extension Ti,r 2 is defined only when Yy and T 2 have disjoint domains. 

In the rest of the paper we assume to work only with well- typed terms. The set of well-typed terms 
will be denoted as T, while P and M denote the set of terms of type sys and mem, respectively. By 
convention we shall use x, y ... for variables of type mem, and X, Y ... for variables of type sys. A 
similar convention is used for base type terms: a, z, ... are terms in M, while P, Q, ... are terms in P. 

Systems can be rearranged according to a structural congruence relation; the intended meaning is 
that two congruent terms actually denote the same system. Structural congruence, =, is the smallest 
equivalence relation on (possibly open) terms of the language which contains the axioms and rules listed 
below. We write T h M = N as shorthand for T h M : t and Fh N :t and M = N. Notice that this notation 
implicitly assumes that structural equivalent terms must be of the same type. 

rh poq = qop rhP°(Q°R) = (poQ)oR rhp°<> = p rhO(ioD = ^ 

rr-a|T = T|a r h <t|(t|p) = (cr]T)]p rr-a|0 = a 
FhP=Q rha=T rhP=2 rha=T 

rhP°R = qor rh o\p = x\p rh gIpd = t<iqd 

rha = /3 Tha = T Tha = T Tha = T 



rha.a = j3.r rh@(a) = @(T) rh »£((t) = ®£(t) 

With respect to the structural congruence of 0, we have added the possibility of rearranging the sub- 
membranes contained in co-phago and pino actions (the last three rules of the table above). 



Reduction Semantics The dynamic behaviour of Brane Calculus is specified by means of a reduction 
semantics, defined by means of a reduction relation ("reaction") ^ C P x P, whose rules are listed in 
Table [2] It is easy to see that subject reduction holds. Note that the presence of (red-equiv) makes this 
not a structural presentation, since term structure can change according with =. 
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; : : : - : : „ — — — (red-pha£o) 

^(p).T\T dQD°^ n .a\a (lPD^T:\T: (lp<la\ao<lPDD QV 

(rcd~cxo^ 

»^.T|T <l8 B .(T|ob<[PD ° 2D — ► ct|<j |t|to<I!2D °P 

@(p).o\o <lPD — ► o\o <lp<L°D °PE> (red " pino) 
P^Q P^Q P = P' P'^Q! Q' = Q 

— „ — — — (red-loc) (red-comp) (red-equiv) 

odPD^odQD PoR^QoR y P^Q 

Table 2: Reduction semantics for the Brane Calculus. 

3 A Labelled Transition Semantics for Brane Calculus 

In this section we introduce a labeled transition system (LTS) for the Brane calculus, along the lines of 
Ifl3l where a LTS is outlined in the case of Mobile Ambients. A SOS style presentation of the semantics 
will help in the definition of the stochastic semantics of Section|4| 

We shall use a meta-syntax for simple syntactic manipulation of terms. The meta-syntax is a simply 
typed A -calculus and can be thought as a primitive system of higher order abstract syntax. We extend 
the base types with function types, in order to type terms in the meta-syntax. We also add to the term 
signature the A -abstraction and application constructs, which should not be considered as a language 
extension: their function is to allow for a structural definition of a labelled transition system, and as such 
they have no computational meaning. 

M ::=... | XX:t.M \ M{M) t ::=... \ t->t 

The A-abstraction binds variables. Terms are taken up-to a-equivalence on bound variables. 

Terms in the meta-language are typed using the standard typing rules for typed A -terms, and the two 
rules below are added to the set of of type rules of Table [T] 

r,X:;hM:/ fhM:Mf' T\~N :t 



rhAX:;.M:^' (lambda) r\-M(N):f ^ 

We consider terms to be syntactically equal up-to j8t] -equivalence, i.e., the smallest congruence 
that contains the two axioms (XX : t. M)(N) = M{N/X} and XX : t. M(X) = M. Thus, for instance, 
(AX : sys. cr<lXD)(P) and adPD are the same term. In the remainder of the paper, when we write a 
meta-syntax term of base type (sys, mem or act), we mean the term obtained by complete j3rj -reduction 
(which exists because our metalanguage is a simply typed-A calculus Q). This assumption will ease the 
description avoiding many technicalities which are out of the scope of the paper. 

We extend structural congruence to meta-syntactic terms by adding the following rule. 

r,X:?hP = <2 



Vr XX :t.P = XX :t.Q 

Structural congruence is compatible with A -terms evaluation in the following sense: 
Lemma 3.1. IfM = N:t^t'andT = S:t then it follows that M(T) = N(S) : t'. 
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(S-pref) 



s„.ct -A a 



(S-pref) 



(S^-pref) 



(» ± -pref) 



adPD ^AAZ.Z(ct'(IPD) 



(phago) 



CT^W 



(phago) 



CT^Act' 



(exo) 



0„(p).ff XT 



(©-pref) 



(7 —XT' 



CT T — > CT T 



(L-par) 



a , 
CT — XT 



p Phag °"; p 



PoQ^i^Xz. (F{Z)°Q) 



U I I 

X CT — > T CT 



(L° -phago) 



(R-par) 



P^A 



PofiJ^jtf.^og) 



(L° -phago) 



P^AS 



o-apD ^A AXy. a'lydXDop 

phago,, phago,, exo„ 
(id-phago-L) 



(L°-exo) 



Act' 



po e i4p (A ) 



id 



PAP' 



CTdPD AS(o)(ct') 

(id-loc) 



(id-exo) 



PAP' 



ct^Act' 



id 



(id-pino) 



ctCPD A CTflP'D 



id 



CTdPD A CT'dpCloDoPD 

— (L°-id) 



P G A P' o g 



Table 3: Labelled transition system (LTS) for Brane Calculus. By convention F : (sys — > sys) — > sys, 
A : sys — >• sys and S : sys — > mem — » sys. Symmetric rules (R°-phago), (R°-phago), (R°-exo), (R°-id), 
and (id-phago-R) are omitted. 



We can now define a labelled transition system for the Brane Calculus, following ITT31 . This system 
will be structural (i.e., in SOS format) and finitary branching; these features will allow to define smoothly 
the stochastic semantics in Section [4] To our knowledge this is the first finitary structural operational 
semantics for the Brane Calculus. (An LTS for Brane Calculus has been given in (H, but it was not 
structural nor finitary branching.) The rules of our LTS are given in Table [3} and are organized into two 
parts: rules for membrane terms and for system terms. In the definition we have omitted types since 
the upper/lower case notation for sys/mem typed variables makes the presentation clearer. The only 
exception to this notation is in (phago) and (L/R°-phago) rules, where the variable Z has type (sys — > sys). 

An example of how labelled transitions are derived using the rules in Table [3j is shown in Figure [TJ 
the derivation leads to the (red-phago) reaction of Table [2| note, in particular, how the (id-phago) rule 
is applied. Similar derivations hold for all rules in Table [2] except for (red-equiv) (which makes the 

reduction semantics "not structural")- Hence, it is easy to see that AC The converse follows by an 
inductive analysis of the structure of the source processes of an /<i-transition. 

Proposition 3.2. IfP A Q then P^Q.IfP^Q then P A 2' for some Q' = Q. 

Labelled transitions are compatible with structural congruence in the following sense: 
Lemma 3.3. IfP A P' and P = Q then there exists Q 1 such that Q' = P' and Q A Q'. 

The labelled transition systems, as usual, induces a bisimulation relation on terms. Due to the gen- 
eralization to meta-syntactic terms, the canonical bisimulation considers also A -terms. Since in Brane 
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(s^-pref) 

« (»-pref) ^±f n \ „ 8 »H „ 

»„.g — » Q ^ (L . par) 

(L-par) » ± (p) 

s„.a|ob ^> a|c7 »,t(P)-^l T o - JL - ^ t|t 

(phago) = (phago) 

v n .o\o Q <lPD ^% AZ. Z(a|a dPD) ^(p).t\tq<IQD ^% AX. T|T dp(LO o gD 

(id-phago-L) 

»„. (y| Ob <IPD ° »^(p).t|t (!2D -h> T|T (!p(!a|ao<lPDD 2D 



Figure 1 : Derivation for the (red-phago) reaction. 



Calculus the dynamics can happen only at the level of system terms P G P, we specialize the definition 
of bisimulation considering only labels in {/J,phago M ,phago„,exo„ | n G A} as follows: 

Definition 3.4 (Strong bisimulation). A bisimulation on Brane Calculus systems is an equivalence 
relation I'CFxP such that for arbitrary P,QeF, (P, Q) G St if and only if 

- ifP P' then 3 Q' such that Q % £/ and (P', Q') G St; 

_ y p phago^ p fhen 3(Jsuch that q i!^so^ q and ^p(g^G(R)) G St for all o,p : mem and R' : sys, 
where R = XX. adpdXD ° R'D; 



ifP A then 3 B such that Q B and (A (R) , B(R) ) G St for all R : sys; 



- ifP T then 3 S such that Q S and (T(R)(p),S(R)(p)) G S for allR : sys and p : mem; 

Two systems P,Q G IP are said bisimilar, written P ~ Q, iff there exists a bisimulation relation St such 
that (P, Q) G St. 

In the definition above, when a transition yields a A-abstracted term M, this term is instantiated 
according to its form. The idea is to recover a bisimulation on only system terms by sufficiently instanti- 
ating the A -abstraction in order to recover a sys-typed term, hence a term in P. This definition is slightly 
different from that in |[T3l for Ambient Calculus, where Rathke and Sobocinski prefer to add to the LTS 
derivation rules for explicit instantiation. Although this design choice is equivalent from the point of 
view of the resulting bisimulation relation, it leads to an infinitely-branching LTS, since they chose to 
endow labels with the instantiation parameters. In the definition of the stochastic semantics for the Brane 
Calculus we need the LTS to be finitely-branching, hence we do not add instantiation rules. 



The next lemma follows directly from Lemma 3.3 



Lemma 3.5. For arbitrary P, Q G P, if P = Q then P ~ Q. 

This ensures that the structural equivalence = is contained in ~, but we can prove also a stronger 
result, that is = Cj ~. Let P = Ofl» ra (IDD; it is easy to see that P ~ <>, however P ^ <>, hence structural 
equivalence does not coincide with strong bisimulation, and in particular ~ equates more terms than =. 



4 A Stochastic Semantics for Brane Calculus 

In this section we present a stochastic semantics for the Brane calculus, following the construction of 
[8]. We assume the reader to be familiar with basic notions from measure theory (see Appendix [A] for a 
brief summary of the used definitions). We start introducing the notation used hereafter and recalling the 
definition of Markov process (MP) and stochastic bisimulation on them (more details are in |f8l). 
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For arbitrary sets A and B, 2 A denotes the powerset of A, and both [A — > B] and B A will be used to 
denote the class of functions from A to B. As usual N, Q, and R denote the sets of natural, rational, and 
real numbers, and Q + and R + the sets of positive rational and real numbers (with zero), respectively. 
Given a measurable space (M,Z), the elements of Z are called measurable sets and M the support-set. 
Let A(M, Z) be the class of measures : Z — >■ R + on (M, Z). From A(M, Z) we distinguish two measures: 
the null measure ft) for which ft) (A) = for all A G Z, and, for a fixed r G R + , the r-Dirac measure on 
N,D{r,N), defined by D(r, iV)(U ie /^f) = Lie/MAf;)* for JV and JV ie/ elements of abase for (M,Z) (thus 
UieiNi G Z), where /n(N') = rifN'=N, otherwise 0. Given two measurable spaces (M , Z) and (AT, 0), 
we use \M — > N\ to denote the class of measurable functions f-.M^-N from (Af , Z) to (A 7 , 0). 

Definition 4.1 (Markov kernels and Markov processes). Le? (M,Z) be an measurable space, and A a 
denumerable set of labels. An A-Markov kernel is a tuple ^ = (M,Z, 8), with 



An A-Markov process of ' M with m G M as initial state, written (^,m), is the tuple (M,Z, 6,m). 

A MP involves a set A of labels which represent all possible interactions with the environment. If 
a G A is a label, m is the current state of the system, and Af is a measurable set of states, the function 
d(a)(m) is a measure on the state space and d(a)(m)(N) G R + represents the rate of an exponentially 
distributed random variable characterizing the duration of the a-transition from m to arbitrary n G N. 

Given a binary relation <% C M x M, we call a subset N CM ^-closed iff & n (N x M ) C N x Af. 
Given (M,Z) a measurable space and ^ C M x M a binary relation over M, with Z(^) we denote the 
set of measurable ^"-closed subsets of M. 

Definition 4.2 (Stochastic bisimulation). For an A-Markov kernel = (M,Z, G) a rate-bisimulation 

relation is an equivalence relation & C M x M such that (m,n) G Si iff for any C G Z(^) an J a G A, 



Two Markov processes (^,m) and (^,n) are stochastic bisimilar, written m n, if m and n are 
related by a rate-bisimulation relation. 

In the rest of the section we define the analytic space of terms, and we show how it can be organized 
as an A-Markov kernel. This will implicitly give a stochastic structural operational semantics such that 
the canonical behavioral equivalence coincides with the bisimulation of MPs. 

Let us define the measurable space of terms. The construction takes place at the level of the meta- 
syntactic terms, hence we will assume to work with terms in T. Let T/= be the set of =-equivalence 
classes on T. For arbitrary M G T, we denote by [M]= the =-equivalence class of M. Note that T/= is a 
denumerable partition of T, hence it is a generator for a a-algebra on T. 

Definition 4.3 (Measurable space of terms). The measurable space of terms (T,Z) is the measurable 
space on T where Z is the o -algebra on T generated by T/=. 

The measurable sets are (possibly denumerable) reunions of =-equivalence classes on T. In the 
following we use T, S, ... to denote arbitrary measurable sets of Z. 

In order to define an A-Markov kernel on (T,Z) we first need to define the set of labels. 

Definition 4.4 (Transition labels). The set of transition labels for ( meta- syntactic) terms is given by the 
pair (A, i), where A = A sys U A mem and 



: A — >■ [M — >■ A(M, Z)J . 



d(a)(m)(C) = d(a)(n)(C). 



(SYSTEM LABELS) 
(MEMBRANE LABELS) 



A sy s — {phago n ,phago n ,exo n | n G A} 
A mem ={»„,S,t(p),»n,Q,t,®n(p) I « G A and p G M} 
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The internal action label id is not in A; we extend A by defining A + = AU {id}. It is no accident 
that the chosen labels are the same of the LTS of Section|3} indeed, the construction of the Markov kernel 
will be guided by the derivation rules listed in Table [3] 

Since we are defining the semantics of a stochastic calculus, we equip the set of actions Act = 
{»„, S^, E>„, E>^ ; @„ | n £ A} with a weight function i : Act — > Q + \ {0}, which assigns to each action 
the rate of an exponentially distributed random variable that characterizes the duration of the transition 
induced by the execution of that particular action. The weight function is such that l(»„) = and 
= this characterizes the fact that two cooperating actions have the same execution rate. 

Now, we aim to define the function 8 : A + —)■[¥—)■ A(T,£)], which will conclude the construction 
of an A + -Markov process for (T,E). To this end, it is useful to give some operations on measurable sets, 
which will ease the exposition of the inductive construction of 6. For arbitrary T,S £ £ and M £ T, 

ns)= U M*r)h i\m= u W\= ?°m= U l N U ^= U l N h 

M-.t^ft'ei N\Mer n°met MdWDet 

N:teS 

The next definition constructs the function 6 : A + — > [T — >■ A(T,E)J by induction on the structure of 
terms. The intuition is that for arbitrary M £ T, 94. £ £ and a £ A + , 8(a)(M)(94) represents the total 
rate of the a actions from M to (elements of) M. 

Definition 4.5. Let G : A + — > [T — > A(T,£)]] be defined by induction on the structure of terms: 
Case 0: For any a £ A + , let 6(a)(0) = CO. 

Case a.M: For any a £ A + let d(a)(a.M) = D(i(act(a)), [M]=) if a = a, d(a)(a.M) = CO otherwise. 
CaseM|N: For any a £ A+ let d(a)(M\N)(1) = 9 (a) (N) (T| M ) + (a) (M) (T^). 
Case XX. M: For any a £ A+ let d(a)(XX. M) = CO. 
Case <>: For any a £ A+ let d(a)(<>) = CO. 

Case M(UVD: For any a £ A mem let 6 (a) (MdNb) = CO. For all other labels in A+ \ A mem 

d(phago n )(M4ND)(?) 
6(phago n )(M4ND)(7) 

6(exo n )(M4ND)(?) 
0(id)(M<LND)(T) 



= 0(»„)(M)([{(7 £ M | AZ : sys^ sys. Z(oWD) £ T}]=) 

= L ^(»»(P))( M )({^ £ M | AX. odpdXD oNb £ T}/=) 
peM 

= 0(»„)(M)({(7 £ M | AXy. a|j(!X!) ° A 7 £ T}/=) 
= B(id)(N)(% m ) + 

neA 

£ 0(©„(p))(M)({(7 £ M | adpfloDoAD £ T}/^) + 
peM 

v 0( g ^)(jV)(5)-e(^)(M)(TQ 

5(KU)(<r' c<r H^«; 
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Case M°N: For any a G A mem let 6(a)(M ° N) = ft). For all other labels in A + \ A mem 

d(phago n )(M°N)(?) = d(phago n )(N)({F : (sys -> sys) sys | AZ. (F(Z) °M) G 7)}/=) + 

6(phago n )(M)({F : (sys -> sys) -»• sys | AZ. (F(Z) ojV) g 7)}/=) 



e(phago n ){M°N){T) = d(phago n )(N)({A : sys -> sys | AX. (A(X) °M) G T)}/=) + 



0(pfeigo n )(M)({A : sys sys | XX. (A(X) °N) G 7)}/=) 
9(exo n )(M°N)(<r) = G(exo n )(N)({T : sys -> mem -> sys | AXj. I(XoM)(y) G 7}/=) + 
0(exo„)(M)({r : sys -> mem -> sys | AXy. r(Xo^)(3;) G T}/=) 
0(id)(M°N)(<r) = 9 (id) (N) (T oM ) + (id) (M) (T oW ) + 

|. A e{phago n ){M){7)-d(ph^3 H ){N){S>L) | 

e(phago n )(N)(T)-dfrha^ n )(M)(A) 

The intuition behind this definition is that each summand corresponds to a derivation rule of the LTS 
of Section [3] For example, in 6(id)(M N)(T) the last summand corresponds to (id-phago-R) rule. 
Similarly, if there are no a-transitions for a term M in the LTS, 6(a)(M) = ft); this is the case of 0, o, 
and XX. M. Note that for each a G M, 6(a) (a) 7^ ft) iff a G act(a). Consequently each infinitary sum 
involved in Definition 14.5 1 has a finite number of non-zero summands. 

In particular we have a correspondence between the LTS and the function 6 in the following sense. 

Proposition 4.6. For arbitrary M £T and a G A+ s the following statements hold 

1. ifd (a)(M) (T) > then there exists M'eT such that M A M', 

2. ifM A M' then there exists M G IT such that M' G T and 6(a) (M) (7) > 0. 

In the proposition above, (1) can be proven by induction on the structure of the term M (assumed 
to be well-typed); while the prove for (2) is by induction on the derivation of M A M' . Note that 
Proposition 4.6 reflects the similarity between Definition 4.5 and the definition of the LTS of Section [3] 



A direct consequence of Proposition |4.6| is the following. 

Corollary 4.7. MAtf iff d(a)(M)([M%) > 0. 

The next theorem states that (T,L,0) is an A + -Markov kernel. Proving this will implicitly show 
the correctness of our construction, indeed it suffices to prove that for each M G T and each a G A + , 
8(a)(P) : £ — > M + is a measure on the measurable space (T,E). 

Theorem 4.8 (Markov kernel for the measurable space of terms). (T,L, 8) is an A+ -Markov kernel. 

A consequence of this theorem is that for each M G T, (T,L, 6,M) is a Markov process, hence we 
can define a stochastic bisimulation for Brane Calculus meta-syntactic terms, simply as the stochastic 
bisimulation of Markov processes in (T,E, 6). 

From the measurable space (T,E) we can recover the measurable subspaces (M,0) and (P,IT), 
respectively for membrane and system terms. Indeed, if is the a-algebra generated from the base 
M/=, and n the one generated from the base P/=, both and IT are contained in E, hence (M,0) and 
(P,n) are two subspaces of (T,E). 
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One may be led to think that by a suitable restriction on 6 we can also recover two A-Markov kernels 
respectively on the sets A mem and A+ s of labels. Of course, it is possible to define two Markov kernels 
by letting mem : A+ em -> [M -> A(M,0)] and s ' ys : A+ ys -> [P -> A(P,IT)J as 

^mem (a) (ff) (*0 = (a) (ff) (iW) for all Stf G 

GsVsO) (*) = (&) ( p ) l 2 ') for a11 * e n. 

Although this definition works well for (M,0), it does not straightforwardly work for (P,IT). In fact 
sys does not enjoy a result similar to Proposition 



4.6 



and moreover the stochastic bisimulation is 
not the equivalence that one may expect. For example, it is easy to see that for all a G A+ s and 
<P G IT, 0(a)(s„d«D)(2>) = 0. The same thing happen for s m doD, that is, 9(a)(»„(I< > D)(fP) = 0, hence 
(P,II, s ' ys , K>„C<>D) and (P,II, s ' ys , » OT d*D) are stochastic bisimilar. By a simple analysis on the structure 
of S n fl<>D and S m d<>D, it i s easy to see that 

e(exo„)(»„doD)([AX;y. 0|jdXD o »)]_) = 0(exo n )(o„(IoD)([AXy. 0|yOX]>)]=) = i(s„) > 
0(phago m )(» m doD)([AZ. Z(0aoD)]^) = 0(phago m )(» m (I^D)([AZ. Z(«)] s ) = i(» m ) > 0. 

We can prove that d(exo„)(v„<loD)('T) > iff XXy. 0\ydXV G T, and 0(phago m )(» m d<>!))(T) > iff 
XZ. Z(o) G T (remember that both 0(exo n )(Q„d«D) and 0(phago m )(» m d«D) are measures). This ex- 
ample points out the problem: the two measurable spaces "are not of the right type". Indeed, in Definin- 
tion 4.5 it easy to see that for labels in A sys (see cases MdND and M °N) 6 contributes to the rate with a 
nonzero value only if the the measurable space given as (last) parameter contains terms of nonbase type. 
In particular, for label phago n the type is (sys — > sys) — > sys, for label phago ;i is sys — > sys, and for label 
exo„ is sys — > mem — > sys. 

This simple consideration induces the definition of another function: 

Definition 4.9. Let d sys : A+ s ->• [P -> A(P,n)] be defined as follows, where PgP and 2> G IT. 

d sys (id)(P){'P) = B{id)(P)('P) 
e sys (phago n ){P){T) = d{phago n )(P)(!F) 

where 7 = {F : (sys — > sys) — > sys | 3 A : sys — > sys. F{A) G 2>} /= 

6 sys (phago n )(P)(P) = d(phago n )(P)(A) 

where M. = {A : sys — > sys | 3F : (sys — > sys) — > sys. F(A) G 2>} /= 

d sys {exo n ){P){?) = 6(exo n )(P)(s) 

where S = {S : sys -)• mem — > sys | 3 <2 G P, a G M. 5(g)(cr) G 2>}/= 

Again (P,L, 6> sys ) is an A+ s -Markov kernel, and in particular the stochastic bisimulation on Markov 
processes is the one that we need. It is easy to see that (P,IT, 6> sys , Q„d*D) and (P,IT, 6> sys , 5£> m floD) are no 
more stochastic bisimilar (assuming ^ l (»)«). 

However, the above definition is not defined inductively on the structure of terms. To this end, we 
provide another function d sys : A+ s — > [P — >• A(P,n)]], defined by induction on the structure of system 
terms, and we prove that d sys is an alternative characterization of 6> sys , i.e., 8 sys and 6 sys coincide. 

Definition 4.10. Let sys : A+ s — > [P — >■ A(P,n)] &e defined on the structure ofP G P, as follows. 

Case P = o; For any a G A+ s , /e? sys (a) = co. 
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Case P = edQb: For any ¥ G II, 

d sys {phago n )(a<lQD)((P) = O mem (v n )(o)({o' 1 rdpda'dQU °RD G !P}/ S ) 
sys p^J(a(!GD)(!P)= £ mem (^(p))(a)({a , |a , (!p(!T(!/?!)!)o ! 2i) G! p}/^) 

peM 

sys (exo„)(CT(ieD)(fP) = e mem (s»)(ff)({<7' | <t'|t(WD ° G 4-}/=) 

e sys (tt/)(ff<[gD)(2>) = Qsys(MO(G)(fP ff n) + 

nsA 

£ emem(@„(p))(a)({a / 1 a' <Ip <[<>]> o <2» G *}/_) + 

peM 

y emem(a,t)(CT)({CT / I g'| ^OgP ° /? € g>}/=) • 6> S ys (eXQ n ) (Q) (?) 

Case P = Q°R: For all a G A sys W 4 G II 

Ma)(0°«)OP) = 0sys(«)(/?)(!Po G ) + sys (a)(e)(!Po R ) 
O sys (id)(Q°R){P) = sys (id)(R)(P oQ ) + sys (id)(Q)(¥o R ) + 



y d sys (phago n ){R)(V) ■ e sys (phago n )(Q)(P) 

k »(»») 



^ OsysiphagOnK Q)^) ■ e sys (phago n ){R){V) 



neA 



The next theorem states that sys precisely characterizes the function 6> sys defined before. 

Theorem 4.11. For a// a G A+ s , P€F,and(P£ U, d sys {a)(P){(P) = sys (a)(P)(fP). 

A consequence of the previous theorem is that (P,n,0 S ys) is an A+ s -Markov kernel, and that the 
stochastic bisimulation on Markov processes coincides with that induced by (P,II, 6 sys ). 

5 Stochastic Structural Operational Semantics and Bisimulation 

In this section we introduce the stochastic structural operational semantics for the Brane Calculus, with 
the aim of defining a behavioral equivalence on system terms that coincides with their bisimulation as 
Markov processes on (P,II, 8 sys ). Notably, it is directly induced from the definition of the function 6 sys 



(Definition 4.10), following the pattern of [8]. Since this is an unusual construction for a structural oper- 
ational semantics, some preliminary discussions are needed. Typically, a structural operational semantics 
for a stochastic process algebra associates a rate r with transitions: P -^> r P' . Instead, in (H transitions 
are not between two processes, but from a process to an infinite measurable set of processes. In order 
to maintain "the spirit" of process algebras Cardelli and Mardare replace the classic rules of the form 
P A P' with rules of the form P — )■ il, where ii is a function from action labels to measures on the 
measurable space of processes. Let us see how this construction can be applied to the Brane Calculus. 

For simplifying the rules of the operational semantics, we first define some operations on the func- 
tions in A(M, ®) A ™ m and A(P, II) A i s , and analyze their properties. We say that a function \i G A(P, n) A ^ s 
has finite support if A+ s \ /I -1 (ft)) is finite or empty (recall that CO is the null measure). 

Definition 5.1. Consider the following constants and operations on A(M, @) Amem defined as follows. 
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• (o mem : A mem — > A(M,0) defined as (O mem (a) = 00, for arbitrary a G A mem ; 

• For arbitrary a G A mem and O ,p G M, let (X a : A mem — > A(M,0) be defined by 

[eUa) = S D ^ £ ^ ^ = £ (for £ G | n G A}) 

I co ifa^e 

[ e ( P )M«) = ( D( ' (eM<ll - ) *' : £ < P ! (^ e { 8 >„|„ e A}; 

[CO ifa^e(p) 

• For arbitrary G A(M,0) Amem with finite support, a G A mem , a, T G M, am/ G 0, fef the 
function jU CT ® T jU'A mem — >• A(M,0) »e defined by 

(li a (DTH')(a)(M)=n(a)(Mr) + n'(a)(M a ) 

Consider the following constants and operations on A(P, 11)^ ancf A(M, 0) Amem defined as follows. 

• ft) sys : A+ s — > A(P, IT) defined as 00 sys (a) = ft), /or arbitrary a G A+ s ; 

• For arbitrary n G A(P,n) As+ y s , ^' G A(M,0) A " K "' with finite support, a G M, P G P, a«<i (P GYllet 
the function jX@p^l' : A+ s — >■ A(P,n) oe defined by: 

(li@ a P li'){phago n ){'P) = i u , (»„)({a' | xdpda'dPn o QD e <P}/=) 
{H@ a P n')(ph^ n )(tP) = £ n'(^(p))({a' I o'dpMQDDoPD G 

PGM 

Ox@?mO(«Wi.)(*) = M'KXK I a'lraeD °P G 
(H@^n')(id)(^P) =n{id)(!P a n) + £ /x'(©„(p))({ff' | a'dpdoD ° e D G -p}/ s ) + 

pGM 

^ /f'(p£)({g' 1 g^ TgPD ° Q G g}/=) • /x(exo„)(g) 



neA 



l(»n) 



Por arbitrary [l , pi' G A(P, n) As+ y s with finite support, a G A sys , P, <2 G P, ant/ fp G IT, /ef the function 
Pp®qP' ': A+ s -> A(P,n) defined by 

(ylp® Q ll')(a)(<P)=yl(a)(<P Q )+ll'(a)(<P P ) 
(Hp® q n')(id) (fP) = M (id) (%) + /*'(«/) + 

„ [}l(j)hago n )(fP)-]X l {phago n ){'E)) + (n(phago n )(p) ■ ji ' {phago n ){T)) 

neA 



Notice that since and ju' have finite support, each infinite sum involved in Definition 5.1 has a finite 
number of nonzero summands. The next two lemmata prove that the definitions of qCDt, p<8>q, and @p, 
for arbitrary a, T G M and P, <2 G P are correct; they also state some basic properties of these operators. 

Lemma 5.2. The following statements hold. 

1. For arbitrary a,T,p G M and jX 1 G A(M,0) Amem with finite support 

(a) n'oCDtV" = H"t®oV', 
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»£} e €{»£,©„} <t->/i' t^^" 

( zel '°) n (P ref ) ; \ — Tj— (pref-arg) . — — — (par) 



( void ) TTTTS TTTr ( loc ) ~^ — 7; 7~~7; 7T ( com P) 



Table 4: Structural operational semantics for Brane Calculus 



(b) (n'oQv At")cr|T®p M'" = M'ctOtIp (ju" t © p jli'"), 
fcj ii , a (Do(O mem = ii'. 

2. For arbitrary P, Q,R G P a«<i il',il",il m G A(P,II) As+ y s with finite support 

(a) ix'p®qH" = h"q®p}i', 

(b) (ju' P 8)Q h")p\q®r H m = h'p®Q\r (h"q®r h'"), 

(c) n'p®<, (0 sys = n'. 

3. w sys @^w mem = w sys . 

Lemma 5.3. The following statements hold. 

1. For arbitrary a, a', T, t' G M a?ici jit', jit" : A mem — >• A(M,0) with finite support 

(a) for e G {»„,»„, fi 6 A}, ff e t implies [e] a = [e] T , 

/or £ G ®n I n S A}, a = T and p = p' imp/)* [e(p)] ff = [e(p')]t, 
(c) a = o' and T = %' imply \1l ' a <S> z \l" = ju'^CDt' M". 

2. For arbitrary P,P',Q,Q' GP,a,TGM, jU,ju',ju": A+ ys -> A(P,n), W v : A mem -> A(M, 0) w/f/i 

support 

(a) P = Q and a = X imply ju@^v = >i@qV, 

(b) P = P' and Q = Q' imply h' p <&q \l" = \i'p®q /J,". 

The rules of the operational semantics are listed in Table[4j The operational semantics associates with 
each membrane a a mapping v G A(M,0) A " ,t,m , and with each system P G P a mapping ii G A(P,n) As V 
For each =-closed set !M G and each label a G A mem , v{a){M) G M + represents the total rate of an 
a-transition of a to some arbitrary element in M; and for each set (P G n and label a' G A+ s , similarly 
ju(a')(fP) G M + represents the total rate of an a'-transition of P to some arbitrary element in T. 

The next lemma guarantees the consistency of the stochastic transition relation — >, and as a conse- 
quence the consistency of the operational semantics. 

Lemma 5.4 (Uniqueness of the measure). For each o G M and P G P: 

1. there exists a unique /J. G A(M, @) Amem such that G — > /J.; moreover, pt has finite support; 

2. there exists a unique jJ. G A(P,n) s y s such that P — > [X; moreover, ji has finite support. 
This operational semantics can be further used to define various pointwise semantics as, e.g. : 

P^Q iff P^pL and p(a)(]Q]=) = r 
Example 5.5. Suppose l(»„) = r, and T,p G M, Q G P; it is easy to see that 
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phciso 2r 

2. »„.a(!P]) o » ;1 .adPD > TdpdaflPDD ° gD ° »„.ct(IPD 

nhciso r 

3. »„.adPD o »„d<>i) — > T(lp(!a(lPi)D ° ° s n fl<>D 

4. s„.ct<IPD o »„doD phag0 "' r ) s n . a dpi) o TdpdoB ° <2D 

5. »„.adPD o s^(p).T(!2D -^4 TdpdadPD!) ° gD 

Note that in (1) the chosen z, p, and Q are not relevant, indeed the "same" transition holds also 
for different choices of them. At first sight (1) seems curious, but the intuition behind it is that it is 
not important where a cell goes once it is phagocytized, but that eventually it could be phagocytized. 
Examples |5 .5 1 (2 - 4) show that the total rate is correctly summed up, while (5) exhibits an internal action 
and how it influences the total rate. 

The next lemma states that operational semantics does not distinguish structurally equivalent terms: 

Lemma 5.6. Stochastic transitions are up-to structural equivalence: 

1. if a =Z and o — > /I then z — > fX, 

2. ifP = Q and P — > /I then Q — > }X. 

Notice that the converse does not hold in general, that is, if for some P, Q e P, P — > pi' and Q — > \l" , 
and \x' = n", this does not imply that P = Q. For example, let P = Od»„dDD and Q = «, then 

P^ (mi = ( %s @!l»Jo)@° 8 , r o« mem ) 
<2^(At2 = « sys ) 

It is trivial to verify that fi\ = n 2 , however P ^ Q. In fact for each a G A sys , }X\ (a)(tP) = (for all 2>), 
since }X\ is of the form [i[ @^&) mem and its result depends only on &) mem . In order to prove il\(id)(p?) = 
it suffices to verify that (£O sys @ <> "[s n ]o)(j<i)(2'o<!p) = 0, but again it is easy because the result depends 
only on co sys , hence jUi = \i 2 - 

Now, we introduce the stochastic bisimulation for the Brane Calculus as the stochastic bisimulation 
on the Markov kernel (P,IT, 8 sys ). We show that systems which are associated with the same function by 
our SOS ar e bis imilar, and that the bisimulation extends the structural equivalence. 



Lemma 



5.4 



shows that the operational semantics induces a function # : P — > A(P,IT) ^ defined by 
0(P) = /x iff P^M 



The next lemma shows that there is a relation between # and the function 6 sys that organizes P as a 
Markov kernel. Indeed, it reflects the similarity between Definition 4.10 and Definition [57TJ 



Lemma 5.7. If (P, n, 6 sys ) is the Markov kernel of system terms and $ : P — > A(P, IT) As y s is the function 
induced by the SOS, then for any P G P, a € A+ s , and T Gil, 

sys (a)(P)(2>) = ^(P)(a)(2>). 

A direct consequence of the previous lemma is that if our SOS assigns to different systems the same 
function, then they are stochastic bisimilar with respect to the bisimulation on Markov processes. 

Corollary 5.8. For arbitrary P, Q G P, ifP — > fx and Q — >■ fl, then P ~(p,n,e sys ) Q- 
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This guarantees that we can safely define the stochastic bisimulation for the Brane Calculus as the 
stochastic bisimulation on (P,n, 6 sys ), as we do in the next definition. 

Definition 5.9 (Stochastic bisimulation on systems). A rate-bisimulation relation on systems in an 
equivalence relation ^CPxP such that for arbitrary P,Q G P with P — > jj, and Q — > jj, ', 



{P,Q)eM iff n(a)(C)=n'(a)(C) for any C G U{M) and any a G 



-sys 



Two systems P,Q G P are stochastic bisimilar, written P ~ Q, iff there exists a rate bisimulation relation 
& such that (P, Q) G £%. 

The next theorem provides a characterization of stochastic bisimulation stating that « is the smallest 
rate-bisimulation relation on P. 

Theorem 5.10. The stochastic bisimulation relation ~ is the smallest equivalence relation on P such 
that for arbitrary P,Q G P with P — > /I and Q — > jll', 

P~Q iff n(a)(C) = n'(a)(C) for any C 6 n(«) and any a 6 A+ s . 

Example 5.11. For arbitrary a G M 

1. Odad^DD^o 

2. ©„(a)|@„(a)(!*D ^ ® n (o).@ n (o)d^ 

3. ^(0)K«oD ~ ^(0).^doD ~ »£(0)<[oD 



Examples 5.1 1 1) and (3) are peculiar to the Brane Calculus; (2) is interesting because in the seman- 
tics without stochastic features the two processes are bisimilar: ©, ! (a)|@„(a)(! <> i) ~ ®, 1 (a).©„(a)(l <> D. 
The next theorem shows that our stochastic bisimulation behaves correctly with respect to structural 



equivalence; it is a direct consequence of Lemma 5.6 and Theorem 5.10 



Theorem 5.12. For arbitrary P, Q G P, if P = Q, then P w Q. 

In addition to to the result of the previous theorem, notice that « is strictly larger then =, indeed 
from Example |5TTTj ;i) we have that, for any a G M, OdaflDD °, however OdadDD ^ o. 



6 Conclusions 

In this paper we have presented a stochastic extension of the Brane Calculus; brane systems are inter- 
preted as Markov processes over the measurable space generated by terms up-to syntactic congruence, 
and where the measures are indexed by the actions of Brane Calculus. For finding the correct actions, 
we have introduced a labelled transition system for Brane Calculus. Finally, we have provided a SOS 
presentation of this stochastic semantics, which is compositional and syntax-driven. 

Stochastic semantics for calculi of biological compartments (but not Brane Calculus) have been given 
in literature; see (Hill for stochastic versions of BioAmbients and lfT5l for a stochastic 7i-calculus with 
polyadic synchronisation. However, these semantics are "pointwise" and not structural, tailored for 
stochastic simulations using Gillepie algorithm. As shown in Section [5} a "pointwise" semantics can be 
readily obtained from the SOS given in this paper. An interesting future work is to investigate how these 
simulation algorithms and techniques can be adapted to our setting. 

There are several directions for further work. First, we would like to prove that the stochastic bisim- 
ilarity (which corresponds to Markov bisimilarity) is a congruence. Then, we want to extend the theory 
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to cover a notion of "approximate behaviour", in order to measure how much two systems are bisimilar; 
this is important in biological contexts, where usually we can compare only with approximate data (e.g. 
coming from experiments). 

We can consider also to add further constructs to the Brane Calculus, like "bind&release" and repli- 

cation. For the latter, we should add rules like Po]P „ p to the LTS of Table 3 on the stochastic side, 

'p— yp 1 n 
these rules would lead to a new case in Definition [43} which we expect to be a nxed point equation to be 

solved in a suitable domain (e.g., complete metric spaces). 

Finally, we would like to apply the present approach to other measurable aspects; in particular, 

geometric (e.g. volumes), physic (e.g. pressure, temperature) and chemical aspects are of great interest 

in the biological domain. 
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A Some measure theory 

Given a set M, a family E of subsets of M is called a c-cilgebra if it contains M and is closed under the 
formation of complements and (infinite) countable unions: 

1. M6l; 

2. A G E implies A c G E, where A c = M \ A; 
3- {A ; } ;GN C I implies UieN A i G s - 

Since MgI and M c = 0, G E, hence E is nonempty by definition. A a-algebra is closed under countable 
set-theoretic operations: is closed under finite unions (A,B G E implies A UB = AUfiU0U0U-- - GE), 
countable intersections (by DeMorgan's law AdB = (A c U B c ) c in its finite and inifite version), and 
countable subtractions (A,B G E implies A\B =AC\B C G E). 

Definition A.l (Measurable Space). Given a set M, and a o -algebra on M, the tuple (M, E) is called a 
measurable space, the elements ofL measurable sets, and M the support-set. 

A set CI C 2 M is a generator for the a-algebra E on M if E is the closure of CI under complement 
and countable union; we write a(Cl) = E and say that E is generated by CI. A generator CI for E is 
a base of E if it has disjoin elements. Note that the a-algebra generated by a CI is also the smallest 
a-algebra containing CI, that is, the intersection of all a-algebras that contain CI. In particular it holds 
that a completely arbitrary intersection of a-algebras is it self a a-algebra. A a-algebra generated by CI, 
denoted by a (CI), is minimal in the sense that if CI C E and E is a a-algebra, then a (CI) C E. Some facts 
about a generators are that, if CI is a a-algebra then obviously a(Cl) = CI; if CI is empty or CI = {0}, or 
CI = {M}, then a(Cl) = {0,M}; if CI C E and E is a a-algebra, then a(Cl) C E. 

A measure on a measurable space (M, E) is a function ii : E — > R + such that 

1. M(0) = 0; 

2. for any disjoint sequence {Ni} ieI C E with / C N, it holds ix(\J ieI Ni) = Lig/A 4 ^')- 

The triple (M,L,jj. ) is called a measure space. A measure space (M,E,/x) is called finite if ju(M) is a 
finite real number; it is called a-finite if M can be decomposed into a countable union of measurable sets 
of finite measure. A set in a measure space has a-finite measure if it is a countable union of sets with 
finite measure. Specifying a measure includes specifying its domain. If ji is a measure on a measurable 
space (M,E) and E' is a a-field contained in E, then the restriction jj,' of ji. to E' is also a measure, and in 
particular a measure on (M', E'), for some M' CM for which E' is a a-algebra on M' . 



A notable measure is the Dirac measure. If CI is a base for (M,E), N £ CI and r G M + , then the function 

/: Cl->R+ 



r-Dirac measure on N. 

Let A(M,E) be the class of measures on (M,E). It can be organized as a measurable space by 
considering the a-algebra generated by the sets Ns >r = {M £ A(M,E) : n(S) > r}, for arbitrary S G E and 
r > (the support-set is A(M,E) and the a-algebra is o (\J SeI/ / . eR + Ns, r )). 

Given two measurable spaces (M,E) and (N,&) a mapping f-.M^-Nis measurable if for any 
T G 0, f~ l (T) G E. Measurable functions are closed under composition: given f-.M^-N and g: N — > O 
measurable functions then g o / : M — > O is also measurable. 




can be extended, by f(\Jiei^i) = LiG//(^')' to a measure on (M,E) denoted by D(r,N) and called the 



